Coexistence of the long-range and short-range magnetic order components in SrEr204 
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Single crystal neutron diffraction reveals two distinct components to the magnetic ordering in 
geometrically frustrated SrEr204. One component is a long-range ordered k = structure which 
appears below Tjv = 0.75 K. Another component is a short-range incommensurate structure which 
manifests itself by the presence of a strong diffuse scattering signal. On cooling from higher tem- 
peratures down to 0.06 K, the partially ordered component does not undergo a pronounced phase 
transition. The magnetic moments in the long-range commensurate and short-range incommensu- 
rate structures are predominantly pointing along the [001] and [100] axes respectively. The unusual 
coexistence of two magnetic structures is probed using both unpolarised and XYZ-polarised neutron 
scattering techniques. The observed diffuse scattering pattern can be satisfactorily reproduced with 
a classical Monte Carlo simulation by using a simple model based on a ladder of triangles. 

PACS numbers: 75.25.-j, 75.47.Lx, 75.50.Ec, 75.40.Mg 
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I. INTRODUCTION 

The geometric frustration of magnetic interactions fre- 
quently results in a highly degenerate ground state of the 
system. In many frustrated systems such as those based 
on the the pyrochlore structure, kagome and garnet lat- 
tices, this degeneracy is often responsible for coopera- 
tive paramagnetism, in which the magnetic correlations 
remain short ranged and the system fluctuates contin- 
uously in the degenerate ground state manifoldii"— The 
degeneracy may be lifted, and a particular type of long- 
range order stabilised at lower temperatures, by thermal 
or quantum fluctuations (through the order-by-disorder 
mechanism). Long-range ordering in such systems is fre- 
quently non-trivial, may be incommensurate with the 
crystal structure and may not prohibit a disordered com- 
ponent remaining down to very low temperatures. The 
coexistence of magnetic ordering and cooperative param- 
agnetism has been recently observed in, for example, the 
disordered system with the fluorite structure TbsNbOyj^ 
It is, however, unusual to observe two different types of 
magnetic order within a single system which does not 
show any signs of chemical instability or phase separa- 
tion. We report such an observation in the rare-earth 
oxide SrEr204, in which a long-range ordered magnetic 
structure with propagation vector k = coexists with a 
shorter-range incommensurate ordering that has a pro- 
nounced low-dimensional character. 

SrEr204 belongs to the family of materials with the 
general formula Srin204, where Ln is Gd, Dy, Ho, Er, 
Tm, and Yb. These oxides crystallise in the calcium fer- 
rite structure^ (space group Pnam), in which the mag- 
netic Ln ions are linked through a network of triangles 
and hexagons^ (see Fig.[T]). When viewed along the [001] 
direction, a honeycomb motif of Ln ions is clearly vis- 
ible, while an edge-sharing triangular motif is apparent 
when viewed along a direction perpendicular to the hon- 



eycomb planes. Although triangular antiferromagnets 
are always subject to magnetic frustration, in a honey- 
comb lattice, where the magnetic ions form a 2D network 
of edge-sharing hexagons, nearest-neighbour interactions 
alone do not result in a frustrated system without compe- 
tition from next-nearest-neighbour interactions. Recent 
theoretical studies^i^ motivated by measurements on the 
antiferromagnetic honeycomb systeni^ Bi3Mn40i2(N03) 
have not reached any firm conclusions as to the nature of 
the ground state in such a system with both nearest and 
next-nearest-neighbour interactions. 

In the absence of reliable data on the strength or the 
sign of the interactions in SrLn204, it is not a priori clear 
whether the correct description of such systems is a 2D 
honeycomb, as quasi-lD chains, or a ladder consisting of 
triangles. A system of similar geometry, /3-CaCr204, has 
recently been shown to behave as a network of triangu- 
lar zig-zag ladders As will be shown below, our Monte 
Garlo simulations suggest that a simple model based on 
a ladder of triangles with a particular ratio of nearest 
and further neighbour exchange interactions captures ad- 
equately the main features of the magnetic diffuse scat- 
tering observed in SrEr204 at low temperatures. 

One important structural feature common to all 
SrLn204 compounds is the presence of two inequivalent 
sites for the magnetic rare-earth ions in a unit cell (de- 
noted as Eri and Er2 for SrEr204 as shown in Fig. [T|). 
The identical ions are coupled to each other through the 
shortest distance, 3.39 A forming chains along the c axis, 
the intra-chain exchange constant is labelled as Ji. The 
chains of identical ions are then coupled into ladders by 
the exchange integrals J2 (Eri-Eri separation is 3.48 A) 
or J3 (Er2-Er2 separation is 3.51 A) shown in red and or- 
ange in Fig. [TJ The ladders are separated from each other 
by the significantly longer distances of 3.87 A and 4.02 A. 
The corresponding bonds are shown in Fig. [T] in grey and 
the exchange interactions are labelled as J4 and J5. The 
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FIG. 1: (Colour online) (a) Positions of magnetic Er"'"^ ions 
within the SrEr2 04 crystal structure. Two inequivalent Er 
sites are shown as dark, Eri, and light, Er2, spheres. The 
chains of identical Er^"*" ions run along the c direction, with 
the shortest distances, Eri-Eri and Er2-Er2, along the chains. 
They are coupled to parallel chains made from either the same 
or different Er'^"'" sites. Different colours used for labelling the 
inter-chain bonds highlight the four distinct links between the 
Er'^^ ions, (b) Ladder of Er2 ions viewed along the [210] di- 
rection (normal to the ladder) showing the nearest-neighbour 
interaction Ji and the next-nearest neighbour interaction J3. 



lattice constants of SrEr204 at 15 K are a = 10.018 A, 
b = 11.852 A, c = 3.384 A (Ref. [111). 

Using powder neutron diffractionfii wc have previ- 
ously shown that SrEr204 orders antiferromagnetically 
at Tjs! = 0.75 K. This transition to an ordered phase 
is also marked by a A-anomaly in the specific heat at 
T/v. The k = magnetic structure below this temper- 
ature was found to consist of ferromagnetic chains run- 
ning along the c axis, with adjacent chains arranged an- 
tiferromagnetically. The refinement suggested that the 
moments point along the c direction, however, only one 
of the two Er'^^" sites possesses a sizeable magnetic mo- 
ment (of 4.5 /Kb at 0.55 K) contributing to this order- 
ing. Within the model constructed from this refinement, 
the magnetic moments may be swapped between the 
two Er^+ sites without changing the calculated powder 
diffraction pattern significantly. 

We have used high-quality single crystal samples for 
the current neutron diffraction study, allowing us to es- 
tablish the presence of another magnetic phase which co- 
exists with the previously reported structure over a wide 
temperature range. 



II. EXPERIMENTAL DETAILS 

Single crystals of SrEr204 were grown by the float- 
ing zone technique using an infra-red image furnace, as 
previously describcdJ^ The three different crystals used 



ranged in mass from 0.5 to 9.5 grams. Wc have per- 
formed neutron diffraction experiments on single crystals 
of SrEr204 using the PRISMA diffractometer at the ISIS 
pulsed neutron source in the UK and the D7 spectrome- 
ter at the ILL in Grenoble, France. The PRISMA instru- 
meni^i^ utilised 16 "^Hc gas detectors with a 1° separation 
providing a good Q resolution, favourable signal to noise 
ratio and high neutron flux in the Q-range from 0.5 to 
5.5 A-i. D7 is a diffuse scattering spectrometer equipped 
with XYZ polarisation analysis^^ which uses 132 '^He de- 
tectors to cover the angular range 4° < 26* < 145° which 
when utilising cold neutrons monochromated to a wave- 
length of 3.1 A results in a Q-space coverage of 0.14 to 
3.91 A~^. These experiments were carried out at tem- 
peratures between 0.05 K and 1.1 K utilising a dilution 
refrigerator inside a standard ^He cryostat. The standard 
data reduction functions were used to create reciprocal 
space maps from the initial data and to normalise the de- 
tector efficiencies from vanadium standards. For polari- 
sation analysis, a quartz standard was used to normalise 
the polarisation efficiency on D7. Background scans were 
performed on D7 using equivalent sample mountings. 



III. NEUTRON DIFFRACTION RESULTS 

The intensity maps shown in Fig. [2] summarise the 
main results of the neutron diffraction experiments. In 
Figs. [2^ and[2lD the intensity in the (hkO) scattering plane 
is shown for two temperatures: 0.06 K, the base temper- 
ature of the dilution refrigerator; and 0.8 K, just above 
T/v. At the higher temperature, broad diffuse scatter- 
ing peaks are clearly present in positions occupied by 
sharp magnetic Bragg peaks below Tjy. Although these 
broad peaks are almost invisible in Fig. [^i when using 
the same intensity scale, they remain present down to 
the base temperature. A highlighted area limited by the 
conditions -1.3 < h < 1.3 and 0.7 < fc < 5.3 in Fig. [2^ 
has intensity ranging from to 0.1 rather than from 
to 1.0 as in the main panels. This area emphasises the 
presence of diffuse scattering (albeit with a significantly 
reduced intensity compared to the Bragg peaks) even at 
the lowest temperature. This weak scattering forms a pe- 
culiar lozenge pattern, which is much more clearly seen at 
higher temperatures (as in Fig. [2b) and which is mostly 
determined by the positions of the magnetic Er"^^ ions in 
the unit cell. 

All of the sharp Bragg peaks observed in the (hkO) 
scattering plane during the single crystal experiments are 
in agreement with the previously reported k = struc- 
ture obtained from the powder refinement. The corre- 
sponding magnetic structure may be found in Fig. 5 of 
Ref. [ill- Example cuts through the (hkO) plane at the 
(120) position are shown in Fig. [3^ for the base tempera- 
ture and 0.8 K, in which the diffuse scattering seen above 
Tn is evident. From the width of the fitted curve corre- 
sponding to the diffuse scattering, the correlation length 
is estimatedi^ as 20 ± 1 A. 



3 




0.6 



0.4 



1.5 



0.5E 



-1 1 2 3 4 5 
(hOO) 


-1 


1 2 3 4 5 
(hOO) 






MiM»»iM» ■mnnfT — 1 — n 






• -Ii 


1 ■ 

• 


il 
■f-i 



-1 0123456 -1 0123456 
(OkO) (OkO) 



0.8 



0.6 



0.4 



10.2 
n 



FIG. 2: (Colour online) Reciprocal space intensity maps of the magnetic scattering from SrEr204 in the (hkO) plane (top panels) 
and in the (Okl) plane (bottom panels) at 0.06 K (left panels) and 0.8 K (right panels). The highlighted area —1.3 < h < 1.3, 
0.7 < k < 5.3 in panel (a) has 10 times lower intensity scale to emphasise the presence at the lowest temperature of a weak 
diffuse scattering otherwise obscured by the much more intense Bragg peaks. The magnetic scattering is isolated from the 
nuclear and spin-incoherent contribution by full XYZ polarisation analysis using D7 diffractometer for the (hkO) plane. In the 
case of maps of the (Okl) plane, the intensity shown is obtained by removing the nuclear contribution from the non-spin-flip 
measurement with neutrons polarised orthogonal to the scattering plane,— following Ref. [l4l . The bounded boxes in all four 
panels delineate the areas of intensity displayed as cuts from the equivalent PRISMA data in Fig. [3] 



Similar measurements taken for the (Okl) seattering 
plane reveal eompletely different diffraction patterns. At 
the higher temperature of 0.8 K, broad diffuse scatter- 
ing intensity forms clearly visible rods, which run along 
the horizontal k axis and which are positioned around 
half-integer points along the vertical I axis as shown in 
Fig. On cooling the system down to 0.06 K, these 
broad diffuse scattering rods became more intense (see 
Fig. [2};). As these rods extend along the k direction, 
the magnetic moments are only weakly correlated along 
the b axis. The short-range magnetic correlations are ei- 
ther two-dimensional (confined mostly to the ac plane) or 
one-dimensional (along the c axis) in nature depending 
on whether the diffuse scattering is limited mostly to the 
(Okl) plane or extends significantly out of this plane. We 
will return to the discussion of this issue after present- 
ing the results of Monte Carlo simulations in the next 
section. 



As can be seen in Fig.[2t, there is an intensity minimum 
around k = 2, however, there is no pronounced decay of 
intensity with increasing k or I significantly more severe 
than that due to the magnetic form factor. The lack of 
marked modulation in intensity of the rods with I and k is 
indicative of magnetic moments directed predominantly 
along the a axis. The positions of the intensity maxima 
along the I axis are clearly incommensurate, however, 
their proximity to the half-integer points suggests that 
the unit cell for this magnetic structure is nearly doubled 
along the c axis. The shape of the rods also changes on 
cooling. The two rods are at positions (0, k,^ + d) and 
(0, fc, I — (5) respectively, where S is dependent upon fc, 
with the sign of S alternating for each increasing integer 
k value, and the magnitude of S decreases upon cooling. 
At the base temperature, the correlation lengthsi£ along 
the [001] direction for the rods, taken from the central 
k = position, are 128(2) A for the (0, k,^ + S) rod 
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FIG. 3: (Colour online) Cuts through reciprocal space, with 
intensity plotted against Q, taken from experiments on the 
PRISMA diffractometer. (a) The (120) peak measured in 
the (hkO) plane at the base temperature of 0.06 K and at 
0.8 K (just above Tjv). The fitted curves shown as solid lines 
are a single sharp Gaussian peak refiecting a strong magnetic 
scattering below Tjv and a combination of a sharp and broad 
Gaussian peaks reflecting nuclear contribution and magnetic 
diffuse scattering above Tjv. (b) The intensity measured along 
the (0, k, 1+5) feature in the (OfcZ) plane at four temperatures. 
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FIG. 4: Temperature dependence of (a) integrated intensity of 
the magnetic (030) Bragg peak as measured on the PRISMA 
diffractometer; (b) integrated intensity of the (0, 0, ^ + S) 
broad feature limited to the region |fc| < i as measured in 
the (Okl) plane on the D7 spectrometer; (c) the value of S. 



and 70(4) A for the (0, fc, | - 5) rod. Only at the lowest 
temperatures, when S is smallest, do the rods begin to 
tend towards the perfectly straight line shape. Fig. [31d 
shows cuts across the {Okl) plane corresponding to / = ^, 
containing one of the diffuse scattering rods, at 0.06 K, 
0.4 K, 0.6 K and 0.8 K. The positions and areas of these 
cuts are indicated in Fig. [5] as bounded boxes. 

Indirect evidence for the existence of low-dimensional 
diffuse scattering at positions of wave vector Q which 
cannot be accounted for in a k = model is also found 
on re-examination of the powder diffraction resultsJ^ 
In this data additional diffuse intensity is observed at 



20 21°and 2Q w 34°, corresponding to the most in- 
tense features seen in Fig. [5}: at the k = positions of 
the two rods. 

The temperature dependence of the intensity of the 
(030) peak, associated with the k = structure, clearly 
demonstrates the appearance of magnetic ordering below 
Tn = 0.75 K, as shown in Fig. S^. Figs. |3Jd and Ht on 
the other hand show that for the diffuse scattering ob- 
served in the (Okl) plane no well-defined phase transition 
temperature is apparent. For example, the integrated in- 
tensity of diffuse scattering measured in the (Okl) plane 
for the k < ^ region at I ^ ^ gradually decreases with 
increasing temperature (see Fig. IJb), while the positions 
and shapes of the rods in reciprocal space, as charac- 
terised by the parameter S, are also smooth functions of 
temperature (see Fig.HJ:). 



IV. MONTE CARLO SIMULATIONS 

In order to gain further insight into the intriguing 
physics of SrEr204, we have simulated the partially or- 
dered magnetic arrangements and associated neutron 
scattering patterns of this compound via the Monte Carlo 
(MC) methodi^i using a standard Metropolis algorithm. 
In addition to the isotropic exchange interactions, a sin- 
gle ion anisotropy term^ was introduced into the simula- 
tion with the b axis being a hard direction of magnetisa- 
tion as suggested by the magnetisation measurements 
The simulations were performed with either 10 x 10 x 30 or 
20 X 20 X 60 supercells containing 12000 or 96000 Heisen- 
berg spins with open boundary conditions>i^ Simulated 
annealing was performed from a temperature T/ Ji = 400 
(to ensure the initial acceptance rate always stayed above 
90%) down to T/Ji = 0.001 with an exponential cooling 
rate of 0.96 and 8000 spin flips per temperature. 

Powder diffraction resultsii showed that only one of 
the two Er^"*" sites undergoes long-range magnetic order- 
ing. One can therefore conjecture that the observed low- 
temperature diffuse scattering is mostly due to the other 
site. In the following simulations we presumed that the 
Eri sites order while the Er2 sites do not. An alternative 
arrangement where the Er2 sites order while the Eri sites 
do not was also examined and gave matching results. In 
the minimalistic model used to simulate the diffuse scat- 
tering, the magnetically ordered Eri sites have been dis- 
regarded and only the Er2 ions and interactions between 
them are considered. The Er2 sites form triangular lad- 
ders running along the c axis that are parameterised with 
two interactions Ji and J3 along the legs and the rungs 
respectively as shown in Fig. [TJ 

Such parameterization produces a network of uncor- 
related ladders where the spins are rotating in the a-c 
plane in each ladder, and the associated scattering func- 
tion displays noticeable modulation around the I = i 
and / = I positions in the (Ofc^) scattering plane in broad 
agreement with the experimental data. The correct mod- 
ulation was obtained for a ratio J1/J3 ~ 5 with both in- 
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FIG. 5: (Colour online) Monte Carlo simulations of the recip- 
rocal space maps of the magnetic scattering in the (Okl) plane 
(top panel) and in the (hkO) plane (bottom panel). The in- 
tensity of magnetic scattering is calculated by averaging 10 
different spin-configurations annealed from T/Ji = 400 down 
to T/Ji = 0.001 (see text for detail). 



teractions made antiferromagnctic in order to match the 
specific shape of the diffuse scattering rods in the (Okl) 
plane. This choice of the interaction signs makes sure 
that the scattering is more intense slightly above (below) 
the I = i position for fc = (fc = ±1), but the situation 
is reversed around ^ = §■ 

To progress further in our analysis, it is useful to note 
the equivalence of the well-studied linear chain model 
with nearest and next-nearest interactions^^ and the Er2 
ladders described here if these were "stretched" along 
the c axis. The observed deviation of the most intense 
diffuse scattering in the {Okl) plane from the (0,0,1/2) 
position translates into a deviation of approximately 6 
degrees from a perfect antiferromagnctic arrangement of 
the nearest spins in the ladder. 

Whilst this simple two-interaction model seems to cap- 
ture most of the physics, some detail, such as the absence 
of scattering around the k = 2 positions both for the 
I = ^ and I = I rods of scattering in the (Okl) plane 
and the very peculiar pattern of intensity in the (hkO) 
plane, are not properly accounted for. In particular, the 
lozenge pattern of weak diffuse scattering, emphasised in 



the highlighted area in Fig. [2^, is completely blurred in 
the simple J3-J1 model with the intensity not confined 
mostly to its corners. This calls for the addition of an- 
other effective interaction between the ladders. 

Keeping J1/J3 ~ 5 constant, several such interactions 
were tried. The best agreement between the simulated 
and experimentally observed patterns of diffuse I{Q) has 
been obtained by an addition of the interaction denoted 
as J' in Fig. [TJ Within the limitations of the minimal- 
istic model used in the simulations, this interaction cou- 
ples the magnetic ions belonging to different ladders and 
should be treated as an effective weak inter-ladder cou- 
pling rather than a specific exchange between ions in cer- 
tain positions. The value of J' cannot be determined 
accurately, as its inclusion influences primarily the over- 
all shape of the diffuse scattering rather than a partic- 
ular feature whose position or intensity could be linked 
directly to the J' /Ji ratio. Much more accurate exper- 
imental results as well as the simulations performed on 
much larger systems (to improve Q-resolution) will be re- 
quired in order to establish J' reliably. The cuts through 
the simulated I{Q) corresponding to the two measured 
scattering planes, (hkO) and (Okl), are presented in Fig. [5] 
for a model in which J' / Ji = 0.01. 



V. DISCUSSION AND CONCLUSIONS 

We return here to the intriguing issue of the dimension- 
ality of the short-range correlations in SrEr204. From an 
experimental point of view, in order to check whether the 
short-range correlations seen as undulated rods of inten- 
sity in Figs.O; and[2Ji are quasi-lD or 2D in nature, one 
has to make extended scans along the a direction across 
the reciprocal space points such as (0, 0, ^+S), (0, 0, |— 5) 
or (0, 4, ■^ + S). These types of scans are impossible on the 
neutron diffractometers used in this investigation with- 
out a complete realignment of the sample. Given the 
neutron beam-time constraints the realignment was not 
performed, therefore the experimental information ob- 
tained so far is insufficient to draw a definite conclusion 
regarding the dimensionality of the correlations. 

When carrying out the Monte Carlo simulations, it is 
straightforward to check whether the intensity observed 
as rods in the (Okl) scattering plane in Fig. [5] protrudes 
significantly out of this plane. Indeed, the results of our 
simulations suggest that this is the case and that the 
intensity forms undulated planes such as (hk^) and (/ifc|) 
rather than rods. One has to bare in mind, however, that 
the simple model adopted in the simulations excludes the 
interactions between the two rare-earth sites. The results 
of the simulations call for further experimental studies, 
which should aim to probe the correlations along the [100] 
direction in order to make a conclusive statement on this 
aspect of the magnetic ordering in SrEr204. 

The onset of the fully ordered k = structure at Tn 
does not appear to affect the diffraction pattern of the 
diffuse phase. In order to explain this lack of coupling be- 
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tween the two types of ordering, and the observation that 
a significant moment contributing to the k = structure 
lies only on one of the two Er'^+ sites, we suggest that 
each of the sites contributes to only one of the two coex- 
isting structures. The different lengths of the inter-site 
exchange pathways (shown in light grey and dark grey in 
Fig. [IJ compared to those between symmetry equivalent 
positions of the same site may provide an explanation 
for the weak correlation between the two magnetic struc- 
tures. 

The low-dimensional character of the magnetic correla- 
tions could be a generic feature of several SrLn204 com- 
pounds. For example, in SrHo2 04 Karunadasa et al^ 
have reported the observation of broad magnetic peaks 
in powder neutron diffraction patterns at low tempera- 
ture. These peaks exhibit a characteristic Warren line- 
shape indicative of the presence of 2D ordering.— It is not 
a priori clear, however, if on further cooling SrHo204 
will undergo a transition into a fully ordered 3D mag- 
netic structure or (similarly to SrEr204) it will main- 
tain separate 3D and low-dimensional magnetic compo- 
nents. In another member of the same family of com- 
pounds, SrDy204, no sharp magnetic Bragg peaks have 
been found down to 50 mK in the recent powder neu- 
tron diffraction experiments;^ while the shapes of the 
observed broad diffuse-scattering peaks are not symmet- 



ric, which again implies 2D correlations. 

In conclusion, we have investigated the magnetic neu- 
tron diffraction patterns in single crystal SrEr204 be- 
tween 0.06 and 1.0 K, and have confirmed the pres- 
ence of a 3D commensurate antiferromagnetic state below 
0.75 K. We have also observed a short-range ordered mag- 
netic phase that coexists with the commensurate phase, 
in which the size of the magnetic unit cell is nearly dou- 
bled along the c axis and which does not exhibit critical 
behaviour. The presence of two non-interacting magnetic 
phases in SrEr204 is an example of the unusual physics 
arising from geometrically frustrated magnetism. 

A classical Monte Carlo simulation shows that a sim- 
ple model based on a ladder of triangles where the near- 
est neighbour interactions arc approximately 5 times 
stronger than the next-nearest neighbour interactions 
mimics satisfactorily the observed diffuse scattering pat- 
terns. Even better agreement between the experimental 
results and simulations is achieved if the ladders are very 
weakly coupled through an effective further-neighbour in- 
teraction. 
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